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Abstract 

Generic low-energy spectral features related to d-d excitations in nearly cubic vanadium per- 
ovskites are predicted for the expected L-edge resonant inelastic X-ray scattering (RIXS) measure- 
ments. Model Hamiltonian describing local electronic properties, including crystal-field effects, of 
vanadium 3d orbitals in the basic t% g configuration is formulated with the help of complementary 
ab initio quantum-chemical cluster calculations. In the presence of 2p-core hole, the local Hamilto- 
nian includes strong 2p-3d electron interactions. As a prerequisite for evaluating RIXS transition 
amplitudes beyond the fast collision approximation, a symmetry-group approach is applied to gen- 
erate a basis set of many-electron wavefunctions of the intermediate core-hole states accessible in 
RIXS processes. Although a comprehensive description of the core- hole multiplets still remains a 
formidable task and requires using specially designed numerical codes, for particular resonant states 
the analysis is simplified and the calculation of RIXS amplitude can be carried out analytically to 
the end. 
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I. INTRODUCTION 



Soft X-ray resonant inelastic scattering (RIXS) is increasingly important in the study of 
electronic properties of strongly correlated 3d systems.- Measurements of d-d excitations that 
correspond to local rearrangements of 3d electrons provide a valuable piece of information 
on the strong coupling between spin, orbital and lattice degrees of freedom in a system. Last 
years, the advances achieved on the instrumental side of RIXS make it possible to measure 
electronic excitations with energy resolution better that O.leV and, hence, to resolve fine 
features of electronic structure not accessible previously. 

The family of vanadium oxides RVO3 (with R=Y or trivalent rare earth ion) with t\ g 
configuration of V 3+ ions have attracted much attention because of its rich and complex 
behavior.- - - In these Mott-Hubbard insulating materials with a weak deviation from the ideal 
cubic perovskite structure the anticipated crystal-field (CF) splitting of t 2g orbitals is rather 
small ~ O.leV and varies with structural changes. In these circumstances, to describe the 
low-energy properties of RVO3 the orbital and spin degrees of freedom have to be treated on 
equal footing, which leads to a spin-orbital superexchange lattice model derived and studied 
in a series of papers." The model was shown to predict diverse behavior of the system 
in dependence on the model parameter chosen and complementary theoretical assumptions 
made. One of the basic assumptions concerns the choice of a point symmetry at V ion, 
the CF parameters and their variation along the RVO3 series irrespective of any orbital 
and spin ordering. We note that CF splitting of t^g orbitals and the multiplet structure 
of t\ g configuration of V ion can be directly inferred from the low-energy part of high- 
resolution RIXS spectra provided the spectral features associated with d-d excitations are 
properly identified as bound excitonic states, i.e., the low-energy in-gap electronic excitations 
produced by RIXS at the final stage in an insulating strongly correlated system. 

In the present study, we examine theoretically the spectrum of d-d excitations with special 
attention to a spin and orbital disordered phase of RVO3. In Sec. II, we begin with a brief 
description of several basic properties of nearly cubic vanadium perovskites, including the 
calculated electronic multiplet structure of the V 3+ ion related to the expected RIXS spectra. 
Fundamentals of the vanadium L-edge (2p-3d core-hole excitation) RIXS transitions are 
discussed in Sec. III. Since a comprehensive analysis of relevant RIXS transitions requires 
knowledge of intermediate core-hole states, this item is considered in Sec. IV in many detail. 
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A generic procedure for calculating the RIXS transition amplitudes together with a more 
detailed consideration of the lowest-energy d-d excitation are present in Sec.V. Summary, 
conclusions, and outlook follow in Sec. VI. 

II. VANADIUM OXIDES WITH NEARLY CUBIC PEROVSKITE STRUCTURE 

At room temperature the crystal lattice of RVO3 is described by an orthorhombic space 
group with the lattice constants a ps b ~ cj a/2. A relatively weak deviation from the ideal 
cubic structure is caused by a regular rotation, tilting and distortion of VC>6 octahedra. As 
a result, the cubic degeneracy of t 2g orbitals occupied by two electrons is removed and the 
accordingly weak CF splitting of ~ O.leV is expected, which is controlled by a variation of 
the R-site ionic radius r R along the RV0 3 series. 

Interplay between weakly split ti g orbitals, spins and lattice degrees of freedom is in the 
hart of a rich behavior observed in RVO3 family, including two types of the orbital-ordering 
(00) and spin-ordering (SO) patterns observed and reported in a number of works, as 
for instance, reviewed shortly in RefU. An ordering, called the C-00, corresponds to the 
orbital arrangement of antiferro-type in ab— plane and ferro-type along the c— axis. The C- 
00 together with the antiferromagnetic G-type spin ordering (G-SO, with spins staggered 
in all three directions) occur at low temperatures, T <100K, in RVO3 with smallest The 
other ordering pattern is the G-00 accompanied with G-SO and observed at intermediate 
temperatures. The G-OO/G-SO is the only ordered phase in RVO3 with larger values of 
tr. The global spin-orbital phase diagram for the RVO3 family including the concomitant 
structural phase transitions derived from the single- crystal measurements was reported.-^. 
Based on a spin-orbital superexchange model for RVO3 a scenario of the puzzling transition 
between the two above-mention ordered phases was suggested.-"- 

In the present paper, our main concern is to the local properties of vanadium valence 3d 
electrons. The local part Hi oc = H c f + H^c of the underlying Hubbard-type model contains 
the crystal field H c f and the Coulomb H d C terms: 

Hd,C = U ^] W m |?l m | + y ' ( U ~Jmm' J n m n m' 2 ^ ] ^mm'S m S m ' + 
m m<m' " m(m' 
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51 Jmm'd mt d ml d m >id m/t . (1) 

Here, is an energy of the m-th 3d orbital, U and J mm i are the Hubbard repulsion and 
Hund's exchange parameters, respectively. When restricted only to the t 2g subspace, the 
Hamiltonian with the same exchange parameters, J m m> — J Hi describes rigorously the mul- 
tiplet structure of t 2g and t\ g configurations.- 1 ^ 

Because the octahedra V06 in RVO3 are tilted and distorted, the exact point symmetry 
at V ion is rather low - the only operation allowed is the inversion, and the cubic 3d basis 
functions, t 2g and e g , are generally mixed. Nevertheless, one expects that a higher point 
symmetry, for instance, D 2 h or D^, can be applied to approximate accurately the CF 
levels and the electronic multiplet structure of the system. In support of this conjecture, 
we have performed complementary ab initio quantum-chemical calculations 11 for a crystal 
fragment containing the rotated and distorted cluster VCV The calculated valence states 
are found to be hybridized V-ion 3d- and oxygen 2p-orbitals, and in a local orthogonal frame 
associated with the principal axes of a rotated octahedron these (antibonding) molecular 
orbitals can be viewed as the canonical t2 g and e g ones. In the following, they are referred 
to as vanadium valence 3d orbitals, and the model Hamiltonian, Eq.(TjQ), is applied to them 
as well. To estimate model parameters, energies of several crystal-field multiplets of the 
valence d 2 configuration, including t 2 and tl g e g , were calculated for the cluster VC>6 and 
compared with those derived from Eq.flT]). We found that the expected for D ih (or D 2h ) 
point symmetry character of a comparatively weak non-cubic CF splitting, both within the 
t 2g and e g subspaces, is well supported by the computed spectra, which is commented in 
detail below. Since the main focus is on the low-energy electronic properties and their 
measurements by RIXS, the calculated multiplet structure of the t 2g configuration is mostly 
discussed. 

To be more precise, we note that in the orthorhombic structure of RVO3 there are four 
non-equivalent crystallographic positions of V-ions coordinated with oxygen octahedra that 
are rotated by the same angles in alternate manner. For a given V ion, assuming the x-, 
y-, and z-axes of the local coordinate system attached to the principal axes of the rotated 
oxygen octahedron, we define the t 2g orbitals as follows: d yz = d%, d zx = d v , and d xy = d$. 
The angular part of the orbital wavefunctions are described by real spherical harmonics: 
Z\ = (i/V2)K + Y?], Z* = (1/V2)K - Y?}, and Z\ = (i/V2)[F 2 2 - Y 2 % with Y 2 being 



4 



the canonical spherical harmonics. 

In our ab initio quantum-chemical cluster approach based on the use of MOLPRO 
computer program^, the wave functions and energies of electron configurations are calcu- 
lated at different levels of accuracy - from a single-reference restricted Hartree-Fock (RHF) 
method through a multi-configuration self-consistent field (MCSCF) ansatz to the multi- 
reference configuration-interaction method (MRCI).— The lattice fragment, whose electrons 
are treated rigorously, contains the cluster VC>6 and eight neighboring R ions, and it is 
imbedded in a large point-charge environment that simulates the Madelung potential on the 
cluster ions. 

In our calculations the primary emphasis was placed upon the most widely studied com- 
pound YVO3. The structural lattice data with the actual rotation, tilting and distortion 
of V06, reported 4 for YVO3 were used. Although the situation is likely to depend on the 
changing degree of octahedral tilting across the RVO3 series^"—, we expect that the generic 
features derived can be directly extended to other members of the RVO3 family, at least 
to those with close values of Tr. Therefore, keeping in mind such an extended usage of 
our cluster calculations, most representative results are reported throughout. As expected 
the calculated subspaces e g and t 2g are well separated by a large cubic CF of ~1.5eV. For 
the high-temperature (T >200K) orthorhombic structure the lowest electronic orbital is d^, 
while the pair of d v , is shifted up by A x ~ O.leV, and these two are only weakly split 
by A2 ~ lOmeV. The inequality Ai ^> A2 can be partly explained for the high-T phase by 
the observation 4 - that in each oxygen octahedron the tetragonal-like contraction of the V-0 
bonds along z-axis is regularly larger than an orthorhombic distortion of the V-0 bonds in 
the x?/-plane. Because the energy resolution of RIXS spectra is usually limited to ~50meV, 
in the multiplet structure fine features due to orthorhombic CF level splitting of A2 ~ lOmeV 
can be hardly resolved. Therefore, in subsequent theoretical analysis of the local electronic 
structure characteristic of the high-T phase of RVO3, instead of D 2 h we use the tetragonal 
D^h symmetry by setting A2 — > and denoting the largest non-cubic CF parameter for the 
t 2g manifold as Ai = A t . This is in contrast to the CF model suggested in RefJo. There, the 
authors assumed the cubic point symmetry Oh with fully degenerate t 2g levels in the high-T 
phase and a symmetry reduction from Oh to D^h (or D 2 h) was associated with structural 
transitions below 200K. 

In high-T phase of RV0 3 the low-energy part of the RIXS spectra is expected to exhibit 
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the multiplet structure of t\ g configuration. Hereafter, a simplifying approximation is used 
that consists in the neglect of comparatively small alternating tilts of oxygen octahedra. 
Thus all the vanadium ions are considered to be equivalent and the 3d orbitals are defined 
in the global (XyZ)-coordinate system with Z-axis along the orthorhombic c-axis, while 
X- and V-axis are rotated by 45° with respect to a and b. 

For the CF of D^h symmetry, the states of t 2g configuration are sorted over the terms 
3 E, 3 A 2 , 1 E, 1 B 2 , 1 Bi, and 2 1 A i , where a subscript indicating the same (even) parity of 
terms is omitted for brevity. For a given term 2S+1 r = (ST), the basis functions \t 2 (ST)M-f) 
labeled by the spin projection M{= —S, —S + 1, S) and the basis index 7 of T, form a 
set of degenerate eigenfunctions of the local Hamiltonian, Eq.(l). The calculated multiplet 
structure is presented schematically in Fig.l. The term energies are found to be: £{ 3 A 2 ) = 
A t , £( l A[) = 2J H - AS', £( l E) = 2J H , £{ l B x ) = £( 1 B 2 ) = A t + 2J H , £{ l A,) = 5J H + A£', 
where A£' m A t [l + (A t / J H )}/3, for A t /J H < 1/2, with the reference energy taken at the 
level of lowest term 3 E of t 2g configuration. From our cluster calculations, it follows also 
that, first, the lowest spin-triplet level coming from the t\ g e l g configuration falls into the same 
energy range - it is stuffed between 1 B\/ 1 B 2 and A%, but not indicated in Fig.l. Second, 
the estimated value of Hund's constant is Jh ~ 0.5eV, which is somewhat smaller than that 
reported usually in literature.— 1 ^ 

For the further purposes, it is instructive to present explicitly the basis functions of 
several terms under consideration. For instance, the ground-state term 3 E comprises six 
states: |t 2 ( 3 £)Mei) = \d vt d (t ), (l/y/2)[\d^d a ) + \d vi d^)], \d vl d a ), for 7 = e 1; and the 
other three states \t 2 ( 3 E)Me 2 ) for 7 = e 2 are obtained by the replacement 77 — > £. The first 
excited term is a triplet \t 2 ( 3 A 2 )M) = Ig^-i-g^), (1 / y/2)[\d^d rj i l ) + Id^d^)}, \d^d v ^), where 
the basis index 7 for the one-dimensional representation A 2 is omitted. In all sets above, 
states are listed according to M = +1,0,-1. With respect to RIXS, the initial (ground) 
and excited (final) states are denoted below as \t 2 , [g]) and \t 2 , [/]), respectively. Here, [g] 
and [/] are the shortenings for the quantum numbers , [g] = ( 3 E)M g e g with M g = 0, ±1 and 
e g = ei >2 , and [f] = (S f T f )M nf . 
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III. RIXS THEORY: BASIC FORMULATION 



Consider a typical scheme of the L-edge RIXS measurement for a strongly correlated 
3<i-electron system whose iV-electron ground state is l^j^)- The incoming photon with 
momentum hk and energy hu^ tuned close to the L-edge absorption promotes an electron 
from the 2p shell to an empty 3c? valence state, thus producing an intermediate core-hole 
state | (\l/ Ar+1 p)[j]). In a subsequent radiative decay of the core hole the emitted (hk', Tiuj-^i) 
photon leaves the 3c? electron system in an excited state with momentum h(k — k') 

and energy hiu^ — uv) that are measured. For given polarization vectors e and e' of the 
incoming and and outgoing photons, respectively, the scattering amplitude of the second- 
order resonant processes governed by the transition operator Dk can be written^ as 

^(k',k;z k ) = ($J ] |(e'D k O t ^(, k )(eD k )|^ ] ). (2) 

Here, Q(z^) is the intermediate-state propagator which describes the system in the presence 
of a core hole: 

y{zk) = 2^ — p , (3) 

where z^ = hu^ + iT, for E\ g j = 0, and the lifetime broadening V of intermediate core-hole 
states is taken to be independent on the state index [i]. We note that T depends on the 
3c?-ion species probed in a measurement and, typically, ranges^ 1 ^ from 0.2eV up to 0.6eV. 
In the dipole limit, the transition operator D k takes the form 

D k = ^E^ kRl E( r i-^), (4) 

where R; denotes the lattice sites of the metal ions and the second summation is over 
electrons belonging to the Z-th ion. The form Eq. (jlj) is especially suitable since we are mainly 
interested in the bound excitonic states, i.e., the low-energy in-gap electronic excitations 
produced by RIXS at the final stage in an insulating strongly correlated system. 
Then, in the second-quantization form the dipole operator Eq.(j4]) reads (Hi = 1): 

D k = 7?E e* kl f £ £<^|^)4A P . + H.c) . (5) 

Here the operator p lmp(T destroys an electron with the orbital projection m p (= 0, ±1) and spin 
a in the 2p shell, while d\ ma creates an extra electron with the same spin in the m-th orbital 



of the 3<i-shell, r = r/|r| , and 1Z 2 contains an integral over radial 2p and 3c? wavef unctions. 
One may check that Dk T = D_k- At this stage, a strong 2p spin-orbit coupling can be taken 
into account by the following replacement 

j,mj 

where Cl™ 3 i are Clebsch-Gordan coefficients for an electron with a total angular momentum 
j— 3/2, 1/2. The insertion of Eq.(jS]) into Eq.© splits the dipole operator: Dk = Dk(j = 
3/2) +D k (j = 1/2). Below we focus on the L 3 -edge scattering processes only and, therefore, 
in Eq.()n]) the sum is bounded to rrij = ±3/2, ±1/2, for j = 3/2. 

By applying the circular components q(= 0,±1) to vectors e and r, the scattering am- 
plitude Eq.pJ) can be rewritten as 

J}';(k',k;z fc )= ^^£^(1/4/;^), (7) 

<?',<? 

where T q > q (e',e) = (— l) q +q e'_ q ,e- q is the polarization tensor and F q t q (g — > f;Zk) is called 
the partial amplitude .— 

From Eqs.d^]),©, and (J7J) one may see that a scattering process being considered as the 
sequence of two local dipole transitions is modulated in space by the phase factor exp[i(k — 
k')l], which makes it possible to measure, besides the local d-d excitation, the dispersion laws 
of accessible propagating electronic excitations. 1 These include spin and orbital waves, and 
their mixture in spin-orbital ordered phases, as well as delocalized excitonic states whose 
propagation may be largely influenced by spin and orbital degrees of freedom. Nevertheless, 
an analysis of the in-gap excitation spectra, together with the corresponding transition 
probabilities and polarization dependence, has to be started with a calculation of a generic 
local dipole transitions. 

The frequently used fast collision approximation^^— may considerably simplify a theo- 
retical analysis of RIXS processes. In general outline, instead of taking proper account of the 
individual scattering events that connect the ground state l^jvi) to a given final state l^m) 
through the intermediary of core-hole states, in this approximation one simply applies an 
averaging procedure that consists in the replacement of the intermediate-state propagator 
by a resonance factor G(zu) — > (<2k ~ ^res) 1 - Such an approximation could be well justi- 
fied if, for instance, the lifetime broadening V of the intermediate core-hole states probed 
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is comparatively large, which is not the case in vanadium oxides where T m 0.3eV is much 
smaller than the intermediate-state interactions. 

To go beyond the fast collision approximation one has to take into account the multiplet 
structure of the intermediate core-hole states explicitly. The problem is extremely compli- 
cated by the fact that in the transition metal ions involved (i.e., V ions) the intra-atomic 
2p-3d interactions are strong, or more precisely, the Coulomb and exchange 2p-3d integrals 
are of the same order as those for the valence 3d orbitals. As known, the multiplet structure 
of core- hole states is directly measured by the X-ray absorption spectroscopy (XAS). Nowa- 
days, XAS spectra for strongly correlated electronic systems are usually predicted and/or the 
measured spectra are interpreted with the use of specially designed computer programs.— - — 

In principle, the computer code for quantum- chemical cluster calculations we used for cal- 
culating the multiplet structure of d 2 configuration of V 3+ ion (Section II), can be extended 
and applied for solving the core-hole multiplets as well. Nevertheless, it is tempting first to 
apply the symmetry-group technique and determine the intermediate-state many-electron 
wavefunctions of interest, and then proceed to the calculation of RIXS amplitudes with- 
out using complementary numerical procedures - the goal that is pursued in the following 
discussion. 

IV. INTERMEDIATE CORE-HOLE STATES IN RVO3 

The first apparent aim is to calculate the on-site matrix elements ( site index 1 is omitted): 
{{d N+1 p)[i]\d\ nu p3 m:j \df g] ) and (df }] \p\ m ,d m > a >\(d N+l 'p) [i] ), which involve two-electron states 

2 2 j 

I^[gj7jyi) °f the low-energy configuration t% g discussed in Section II. The intermediate states 
\{d N+l p)\{\) containing one more electron in 3c? valence shell and a core hole in 2p shell have 
to be determined as eigenvectors of the local Hamiltonian Hi oc extended by including the 
2p-3d interactions H pd ^x and the spin-orbit coupling H PtSO of 2p electrons 

Hioc = H c f + Hd,c + H p d t cx + H PjS0 . (8) 

The strength of the direct Coulomb (C) and the exchange (X) interactions entering the 
Hamiltonian H pdC x are parametrized with the Slater integrals F 2 d , and GL, Gp d , 
respectively. For the 3d 3 2p 5 = (d 3 p) configuration of V ion, representative values^ are: 
F 2 d w 6eV, Gh d ~ 4.4eV and GL ~ 2.5eV. The monopole integral F® d provides a common 
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shift of all the multiplets under consideration. 

To consider the full multiplet structure of the intermediate core-hole states, an appropri- 
ate basis set is build up in few steps. First, we define the auxiliary states |d 3 (S'ir 1 )M 1 7 1 ) of 
the d 3 configuration not disturbed by Hp^cx- Similar to the d 2 configuration, this is done 
by specifying a set of terms (S'iTi), each composed of degenerate eigenstates of H c f + H d C 
with the corresponding eigenvalues Sd(SiTi). At this step, the approximation can be used 
that consists in the neglect in H c f of a weak tetragonal distortion A ~ O.leV compared 
to a strong cubic CF component ~ 1.5eV. Regarding the intermediate states available in 
RIXS processes, this approximation is justified as follows: Like in X-ray absorption, any CF 
distortion, as well as a many-body interaction, that causes small energy effects compared to 
the lifetime broadening V ( w 0.3eV for vanadates) cannot be resolved. Loosely speaking, 
the fast collision approximation can be safely applied separately to each narrow group of 
(quasi-)degenerate intermediate states, but not to the total set of states. 

In case of cubic CF symmetry, the t 2g configuration is given by four (SiTi) terms, namely, 
A A 2 , 2 E, 2 Ti, and 2 T 2 , with the lowest A 2 and degenerate 2 E and 2 T\. With the reference to 
the lowest level £d( 4 A2), the term positions are £d( 2 E) = £d( 2 Ti) = 3Jh and S^T^) = 5Jh- 
For the further purposes, four spin states of 4 A 2 deserve to be presented explicitly (the sole 
7 is omitted): 

\t 3 ( 4 A 2 ),M 1 = 3/2) = |%^ t d ct ), 

\t 3 ( 4 A 2 ), Mi = 1/2) = —(\dtfdrtda) + \dtfd vi d (t ) + (9) 

\t 3 ( 4 A 2 ),M 1 = -1/2) = -^=(\d a d vi dtf) + \d a d^d a ) + \dtfd vi d a )), 
\t 3 ( 4 A 2 ),M 1 = -3/2) = \d a d ril d a ). 

The states contained in other terms (SiTx) can be presented in a similar manner as well. 

Vector coupling of two irreducible representations Si and Ti leads to an equivalent de- 
scription of the t 2g configuration in terms of representations F with the corresponding states 

|t 8 ,P7(5 1 r 1 ))= £ [fis^MaiW^ , (10) 

A/1,71 

where coefficients of the unitary transformation, Eq. (fTUI) . are tabulated 3 ^ (both for 

the single and double- valued point groups). Next we note that for the 2p core-hole with total 
angular momentum j = 3/2, the states \Pa m ) transform according to the four- dimensional 
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double-valued representation G 3 / 2 (T 8 in the Bethe system) of the cubic group, and the 
notation \p 3 ) = \p , G 3 / 2 m,), with m,- = ±3/2, ±1/2, is more convenient. 

— i m i — 

Finally, the complete basis set \t 3 p , Tj'jj) for the core-hole configuration t 3 p is determined 
by vector coupling of the representations T for the "undisturbed" valence t\ g configuration, 
Eq. ( TTUjl . and G 3 / 2 for the 2p core hole. The coupling relation reads 

\t 3 p ,r j7j (rG 3/2 )(s 1 r 1 )) = £ \t 3 7 T^(s l v 1 ))\p ,G 3/2 mj)u^f J , (11) 

7 ,TOj 



re? r 

with the coupling coefficients U^ m 3/ ^ 3 defined as in Refj30l In general, when the pro- 



cedure outlined above is applied, several repeating representations Tj,T'j, T", etc., of 
the same dimensionality may be generated. Repeating representations are discriminated 
with additional labeling, (SiTi) and (rG 3 / 2 ), kept in Eq. lflTI) . to indicate different routes 
(.SiTx) — > (TG 3 / 2 ) — > Tj the representations T j or T'j, Y"j arise. Below we refer to (SiTi) as 
the generating multiplet. 

Eq.t flTT) suggests a conventional basis for a diagonalization of the local Hamiltonian, 
Eq.flE]). Insofar as states of repeating representations are mixed by H p d,cx, the matrix of 
H[ oc takes a block-diagonal form, which can be solved numerically to give the full multiplet 
structure of the core-hole configuration t 3 p. Here, the main difficulty stems from the fact 
that complicated vector-coupled functions, Eqs. (|T0|) and (TlTl) . come into play and calculation 
of numerous off-diagonal matrix elements of Hi oc becomes very tedious. Nevertheless, close 
analysis shows that in several cases calculations are greatly simplified and some prominent 
features of the multiplet structure of a core-hole configuration can be evaluated without 
dealing with the full form of high-rank matrices of H\ oc . A particular example is considered 
in the next Section. 



V. RIXS TRANSITIONS IN RV0 3 

Let us first apply the outlined procedure for examining the states of those Tj's that 
are related to the lowest generating term (Siri)= 4 A 2 . Because A 2 the is one- dimensional 
representation, the Eq.f llOp allows the following identification 

|t 3 ,r"7( 4 A 2 )> = |t 3 ,r = G 3/2 ,7 = M 1 ( 4 A 2 )) = \t 3 ^A 2 ),M x ), (12) 

with the final kets defined in Eq.([9]). Then, according to Eq. dlip . all irreducible representa- 
tions Tj's generated from the term 4 A 2 are contained in the direct group product G 3 / 2 (g)G ? 3 / 2 
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and given^ - in the following decomposition: G 3 / 2 <E> G3/2 = A\ + A 2 + E + T\ + Tj' + T 2 + T^. 
Without 2p-3(i interactions the terms are degenerate, however, being subjected to H pd) cx 
they are split into two subsets in such a way that each subset is coupled by H pdC x to states 
of other multiplets in different manners. To exhibit this splitting, the term energies are 
preliminary found to the first order in H p ^cx by calculating the corresponding diagonal 
matrix elements of the local Hamiltonian, Eq.flH]). This yields three lowest terms, A 2 , T\ and 
T 2 , residing at the same energy that is taken below as the reference (zero) level. Then, the 
other four terms are shifted to higher energies as follows: A\ at j^G pd + ^G 3 d , E and T' 2 at 

4 r 1 J- 6 r" 3 nnrl T> at 17 r 1 4- 2 

15 ^pd 35 ^pdi dllu J l dL 30 7 (j pd- 

The states of three lowest terms, A 2 , T\ and T 2 , the first subset, can be regarded as "edge" 
resonant states for several reasons. Firstly, to any order in H pd> cx the "edge" states are not 
mixed with other states of the t 3 p configuration - the corresponding off-diagonal matrix 
elements of H [oc are checked to be zero. To the contrary, the second subset, i.e., the states of 
complementary multiplets Ax, E, T 2 and T{, are mixed up with states \t 3 p ,Tj n fj(SxTx)) of 
similar multiplets generated from (SiTi) 7^ A 2 and located well above the reference level. In 
total, this mixture gives rise to a broad and dense manifold of intermediate states \t 3 p , Tj'jj), 
but still the energy of "edge" states is retained below the bottom of this broad manifold. The 
"edge" states are stabilized mainly due to, first, that the initial energy of the generating term 
4 v4 2 is well below the positions of other terms (SiTi) ^ A A 2 and, second, for Tj = A 2 , 7\ and 
T 2 the exchange coupling of the valence states |t 3 , Mi^A^) to the core states \p, G^/ 2 rrij) due 
to H P d t cx is optimized. Experimentally, concerning the L 3 -edge XAS/RIXS measurements 
in RVO3, the "edge" states should be associated with the lowest peak expected at the very 
edge of the XAS absorption spectra, while in RIXS they are accessible by tuning the incident 
photons to the lowest resonance energy f)^, i.e., as hco^ — > fifes- 

Thus emerging multiplet structure of the core-hole configuration t 3 p suggests that close to 
the "edge"-states associated with the resonance at Q$ s , the remaining terms Tj evolve into 
a peak-shaped structure of several features at higher energies, f^es > ^res-> f° r ^ = 2,3, .... 
Particular peaks initially assigned to different Tj's are broadened due the core-hole lifetime 
effects. A group of peaks may overlap if their energy spacing is less than the core-hole 
lifetime broadening. This means that a set of core-hole states \t 3 p , Tj7j) that are included 
in a group of "closely spaced" T/s, are not resolved and should be attributed to a peculiar 
resonance energy which is denoted below as {Tj7j}g (n). Then the partial amplitude 
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of corresponding local transitions can be expressed in the following form 

p ( ^fo(n)\ (t 2 , [f]\v qt \t 3 P ,Tju)(t 3 P ,r j7j |P g |t 2 , 

{rj7j}e(n) Zk — nilr 



where V q is the q— th component of the dipole operator, Eq.([H$, taken at k = 0. Eq. (fl3"|) 
shows that the states involved in {Tj7j} e (n) interfere to yield a joint transition at the 
resonance energy hVt^ s . 

To predict the relative positions of f^es; the eigenvalue problem should be solved by di- 
agonalizing numerically the local Hamiltonian, Eq.flH]), with the basis set defined in Eq. flTT]) . 
Together with the calculated eigenstates \t 3 p ; Tj'jj), this would give one a firm ground for 
a comprehensive account of RIXS transition probabilities. As mentioned above, instead of 
performing this ambitious program, here we focus on solving a restricted problem - cal- 
culation of the RIXS transitions that occur via the resonant "edge" states, which can be 
accomplished analytically to the end. 

The explicit form of "edge" states follows directly from Eqs. ffTO]) and ffTTl) . In Eq. ffTTj) 
the coupling coefficients U^ 3 /^ J required for T = G3/2 and Tj = A 2 , T\ and T 2 , are found 



from tables of Refj30|. For the sake of brevity, we reduce slightly the notation as follows: 
^Mim^yj 2 ' 1 " -7 = UmJ^. Complementary reduction includes both |g, G 3 / 2 mj) = \Pi m j)i where 
rrij = ±3/2, ±1/2 for the core-hole states, and \t 3 p; Tjjj(G 3 /2G 3 /2)( 4 A 2 )) = \t 3 p;Tjjj) for 
the "edge" -state kets. In shortened notations the latter can be now expressed as 

l* 3 P,r j7 J>= E |t^(%)M 1 )b,m,K ^ / f4, (14) 

M\,mj 

where the kets \t 3 , ( 4 t4 2 )Mi) are defined in Eq.flH]). 

As — > fl^Jg, the local RIXS transitions are described by partial amplitude, Eq.(Ti"3"l). 
where {Tj7j}g (1) denotes the lowest energy degenerate multiplets Tj = A 2 , T x and T 2 of 
the t 3 p configuration. Matrix elements of the dipole operator are given by 



(t 2 , [f]\v q ,\t 3 p ,r j7j ) = n 2 £ E (t 2 , [f]\d rm7 \t 3 ,M 1 ) 

Ml ,<r m,m p 



■P2 L 



(Y mp \r - q ,\Z m ) , 
(15) 

The complex conjugate of Eq. ffTS"]) and the replacement q' —> q, [f] — > [g] then yields 
(t 3 p , Y j^j\V q \t 2 ^ [g]). Having the states \t 2 , [f]/[g]) and |t 3 ,M!) as defined in Section II 
and Eq.Q, together with the symmetry coupling U^/J^. and Clebsch-Gordan C 3 ^™^ coef- 
ficients known, the transition amplitudes at hu^ — > ttfj are calculated straightforwardly. 
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A connection of the ground state \t 2 , [g]) to particular final states \t 2 , [/]) depicted in Fig.l 
and a polarization dependence of the transitions are of the main interest. 

The calculations show that besides the elastic transition, \t 2 , ( 3 E)M g e g ) — > \t 2 , ( 3 E)M' g e g ) , 
the only excited term of t\ g configuration accessible via the resonance at f2^ s is 3 A 2 , Fig.l. 
This occurs in the processes \t 2 , ( 3 E)M g e g ) ->• \t 2 , ( 3 A 2 )M f ) with the initial and final states 
being degenerate in M g {= 0, ±1), e g {= ei,e 2 ) and M/(= 0,±1), respectively, for the disor- 
dered high-T phase. As seen from Fig.l, the corresponding excitation energy is huf g = A t , 
and thus the basic CF parameter can be measured directly together with a fully identified 
excitonic state. A possible band broadening of this excitonic state is discussed later on. 
Identification of this peculiar d-d transition in RIXS measurements of RVO3 would enhance 
a better understanding of the system. For instance, a trend of interest that A f develops upon 
replacement of R ion in RVO3 can be followed. Second, from Eqs. f|T3|) .f lT5|) one obtains the 
transition intensity 

r'^-A^wj = ay, E EI^(^ e )^(( 3 ^) M ^^( 3 ^)M /; fi«)| 2 ,(i6) 

M f M g ,e g q',q 

where A is the known constant.- Note that in Eq. (fl6|) an averaging over the initial states and 
summation over the final states are performed assuming the state degeneracy mentioned. 
From this, the calculated polarization dependence is given up to a factor as 

r> /s ^A t ,fi«) ~ (4 2 +e ;)e 2 z , (17) 

which may serve for a verification of the previously discussed orbital occupancy both in the 
initial and the excited states. 

The RIXS d-d transitions to the states of higher-energy terms (SfTf) y^ 3 A 2 listed in 
Fig.l can be measured via the resonance processes that occur as hu-k with n > 1. 

Contrary to the lowest resonance, n = 1, now the process \t 2 , ( 3 i?)M 9 ec,) —> \t 2 , (SfT /)M/7/) 
branches off, meaning that for given n(> 1) there are several accessible final terms (SfTf). 
Theoretically, the intertwined transitions involved could be unraveled provided the relative 
energy positions of the resonance energies fl^ and the explicit form of intermediate core- hole 
states Tj7j e (n) are found, which would again require using computer codes. Although the 
overall inspection of the RIXS processes is not given herein, nevertheless a valuable piece of 
information regarding the expected RIXS spectra is predicted, as presented in Fig.l together 
with complementary estimates in Section II, and the lowest d-d excitation is studied in great 
detail. 
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So far, when classifying the multiplets of the t 2 configurations based on the symmetry- 
group technique, we assumed that in high-T phase a deviation from the symmetry is 
too weak to be resolved experimentally With decreasing temperature, one expects that the 
deviation increases and the experimental resolution permits the main effects of the symmetry 
reduction D ih — > D 2h to be observed in RIXS spectra. Although, a reduction of the former 
multiplet structure of t\ g configuration to a new one due to the local symmetry lowering is 
a routine theoretical procedure, a precise analysis of the expected fine structure of RIXS 
spectra at low temperatures is further complicated by the spin-orbital ordering. Actually, 
since in the lattice the neighboring V ions are now coupled by the inter-site spin-orbital 
superexchange, the states of a local multiplet are split by internal exchange fields and the 
former degeneracy is removed. In particular, the symmetry breaking of the ground-state 
multiplet 3 E gives rise to a number of low-energy spin/orbital excitations whose energy 
position (and dispersion properties) are set by the superexchange constant J ~10meV, and 
hence the excitations are expected^^— well below O.leV. From the perspective of RIXS 
measurements, the energy region of interest, huf g <0.1eV, is typically submerged by the 
tail of the huge elastic peak, which makes the selection of the lowest spectral features a 
problematic procedure. 

In this respect, the promising dispersive quasiparticles for the RIXS measurement orig- 
inate from the d-d excitations thus far regarded to be localized excitons, i.e., the bound 
electron-hole pairs in the valence t 2g subshell of V ion. A mechanism behind the exciton 
derealization is very similar to that proposed for the orbital wave propagation.- It involves 
a simultaneous exchange of two active t 2g orbitals at neighboring V-V pairs forming three 
distinguished bonds along lattice axes. For instance, since in the lowest d-d excitation of 
3 A 2 symmetry the d^ orbital involved is inactive along z||c axis, the dispersion law of the 
corresponding exciton is of pure two-dimensional form: -E'(q) = E + (W/4:)(cosq x + cosq y ). 
Here, rough estimate for the bandwidth is W = aW, where W ~ 50meV and the parameter 
ct(< 1) takes into account the temperature dependent spin and orbital intersite correlations 
in different phases. The excitonic gap Eq is determined by the t 2g orbital CF splitting A t 
renormalized compared to A t due to intersite superexchange coupling. In disordered high-T 
phase of RVO3, one has E = A t , assuming that a small self-energy correction due to virtual 
intersite tS^-electron hoppings is already included into A t that should be measured. We 
conclude this short discussion with a general remark that much more work is needed to dis- 
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criminate all possible delocalization mechanisms responsible for a formation of propagating 
excitonic states. It seems that a further improvement of the experimental RIXS resolution 
towards the 10 meV level is required to measure these quasiparticles. 



VI. CONCLUSION 



To summarize, we developed the RIXS theory and applied it to predict most prominent 
features of the low-energy d-d excitations that can be measured in RVO3 - a prototypical 
Mott-Hubbard insulator with t| configuration of V 3+ ions. Complementary ab initio cluster 
calculations allowed us to estimate the crystal field parameters and choose a proper basis set 
for vanadium 3c? orbitals. Initial (ground), intermediate, and final states accessible in RIXS 
processes were classified with the use of the symmetry-group technique, which is prerequi- 
site to the subsequent analysis developed. To go beyond the fast collision approximation 
frequently used for calculations of the RIXS transition amplitude, we examined the inter- 
mediate core-hole states, including the spectral and symmetry properties of many-electron 
wavef unctions, with more detail than usually done before. Eventually, the problem of deduc- 
ing a detailed overall picture for the intermediate states was reduced to a diagonalization of 
a local Hamiltonian which involves the valence electron d-d and the core-hole-valence p — d 
interactions. Mostly because of the huge number of many-body core-hole states entangled, 
the diagonalization procedure is far from being a trivial one. As well known, such a dif- 
ficulty one encounters in any theoretical analysis of XAS spectra^ in strongly correlated 
electron systems, which requires using specially designed computer codes. At this stage, 
we avoided doing any numerical calculations, but instead pursued an analytical treatment 
so far as possible and found, for instance, that a comprehensive analysis of the lowest d-d 
excitations is attainable in a simple fashion. For the basic t\ g configuration, the relative 
spectral position of d-d excitations of higher energy is predicted and expressed in terms of a 
model Hamiltonian describing the local electronic properties of RVO3. In the present study 
we restrict ourselves to the calculation of RIXS spectra expected in a spin- and orbital- 
disordered high-T phase of RVO3. An extension to ordered phases, which is only shortly 
discussed, may be a subject of a further work. Delocalization of low-energy d-d excitations 
leading to a formation of propagating excitonic states was also discussed and some predic- 
tions were formulated. Along this way, a future careful theoretical study and subsequent 
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examination by RIXS of low-energy excitations in RVO3 will provide a better understanding 
of the microscopic mechanisms responsible for rich properties observed in this promising 
system. 

To conclude, we believe that the suggested theoretical approach may be extended and 
applied to a description of selective RIXS transitions in other transition metal oxides with 
partially filled ti g orbitals. In further perspective, the ab initio quantum-chemical cluster 
approach can be also applied first to calculate multiplet structure of both the ground-state 
t^g and the core-hole t^ +1 p configurations and, second, to evaluate the matrix elements of 
the dipole transition operator between them. In this context, a symmetry group analysis 
developed here for a particular case of V 3+ can be immediately developed to provide a com- 
plementary tool for classifying the calculated many-electron wavefunctions and controlling 
the symmetry selections of different d-d transitions. 
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FIG. 1: Multiplet structure of t^g configuration of V 3+ ion in crystal field of D^h symmetry. The 
term branching with the symmetry reduction Oh — > is shown for clarity. In RIXS the initial 
(ground) state is given by 3 E and other terms correspond to local d — d excitations. More precise 
estimates for the excitation energies are given in the text. 
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